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Recent results from jet production in deep inelastic ep scattering to investigate parton 
dynamics at low x are reviewed. The results on jet production in deep inelastic scattering 
and photoproduction used to test perturbative QCD are discussed and the values of 
| a s (Mz) extracted from a QCD analysis of the data are presented. 
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o 1 Introduction 

^ ■ 

Jet production in neutral-current (NC) deep inelastic ep scattering (DIS) and photoproduction 
provide tests of perturbative QCD (pQCD) calculations and of the parametrisations of the 
V ' parton distribution functions (PDFs) in the proton. Jet cross sections allow the determination 
of the fundamental parameter of QCD, the strong coupling constant a s , and help to constrain 
^ ! the proton PDFs. 



Up to leading order (LO) in a s , jet production in NC DIS proceeds via the quark-parton model 
(QPM) (Vq -> q, where V = 7* or Z°, Fig. Et), boson-gluon fusion (BGF) (Vg -> qq, Fig. Eh) 
and QCD-Compton (QCDC) (Vq — > qg, Fig. Ht) processes. The jet production cross section is 
given in pQCD by the convolution of the proton PDFs and the subprocess cross section, 



da ict = dx f a (x, jj Fp ) da a (x,a s (jj R ), fx R , fx Fp ), 

a=g,q,g 

where x is the fraction of the proton's momentum taken by the interacting parton, f a are the 
proton PDFs, fi Fp is the proton factorisation scale, a a is the subprocess cross section and hr is 
the renormalisation scale. 

The main source of jets at HERA is hard scattering in photon-proton (photoproduction) inter- 
actions in which a quasi-real photon (Q 2 w 0, where Q 2 is the virtuality of the photon) emitted 
by the electron beam interacts with a parton from the proton to produce two jets in the final 
state. In LO QCD, there are two processes which contribute to the jet photoproduction cross 
section: the resolved process (Fig. |TJi), in which the photon interacts through its partonic con- 
tent, and the direct process (Fig. [T^), in which the photon interacts as a point-like particle. 

*Talk given at the "XXXIV International Symposium on Multiparticle Dynamics", Sonoma County, Cali- 
fornia, USA, July 26 th - August 1 st , 2004. 
t Ramon y Cajal Fellow. 
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The cross section for the process ep — > e + jet + jet + X is given by the convolution of the 
flux of photons in the electron, the parton densities in the proton, the parton densities in the 
photon and the subprocess cross section: 



da ep ^ ct jet = ^2 / d yfi/e(v) / dx i fih&nPF-i) I dx v fj/p( x P^F p ) da i{l)j (i(>y)j -> jet jet), 
y Jo Jo Jo 

where f 1 / e is the flux of photons in the electron, usually estimated using the Weizsacker- 
Williams approximation, y is the inelasticity variable, fj/ p are the parton densities in the proton, 
determined from global fits, x p is the proton momentum taken by the interacting parton, y,p is 
the proton factorisation scale, /j/ 7 are parton densities in the photon, determined from fits to 
deep inelastic ey data, x 7 is the photon momentum taken by the interacting parton, \ip. is the 
photon factorisation scale, and da^j^i^j — > jet jet) is the subprocess cross section, which is 
calculable in pQCD at any order. 
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Figure 1: Examples of Feynman diagrams for deep inelastic ep scattering processes: (a) quark- 
parton model, (b) boson-gluon fusion and (c) QCD Compton. Examples for photoproduction 
processes: (d) resolved and (e) direct. 

All the data accumulated from HERA Pfl |2] and fixed-target [S] experiments have allowed a 
good determination of the proton PDFs over a large phase space and the evolution of the PDFs 
with the scale hf v for large scales has been successfully described by the DGLAP evolution 
equations jlj. Measurements of jet production in NC DIS jSJ E] and photoproduction |H] 
have provided accurate tests of pQCD and a determination of the fundamental parameter of 
the theory, a s . Most of these measurements refer to the production of jets irrespective of their 
partonic origin -quarks or gluons- and, therefore, have provided general tests of the partonic 
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structure of the short- distance process and of combinations of the proton and/or photon PDFs. 
The identification of quark- and gluon-initiated jets would allow more stringent tests of the 
QCD predictions. 

At high scales, calculations using the DGLAP evolution equations have been found to give a 
good description of the data up to next-to- leading order (NLO). Therefore, by fitting the data 
with these calculations, it has been possible to extract accurate values of a s and the gluon 
density of the proton. However, for scales of Ej, ~ Q, where E^ is the jet transverse energy, 
large values of the jet pseudorapidity, r]^ et , and low values of x discrepancies between the data 
and the NLO calculations have been observed. This could indicate a breakdown of the DGLAP 
evolution and the onset of BFKL |5] effects. These discrepancies could also be explained by 
assigning a partonic structure to the exchanged virtual photon or a large contribution of higher- 
order effects at low Q 2 . 

This report includes recent results of the studies on parton evolution at low x from the HI 
Collaboration, namely forward-jet cross sections as a function of x and the azimuthal correlation 
between the hard jets in dijet events in NC DIS, and the tests of pQCD at high scales from the 
ZEUS Collaboration, namely dijet and three-jet cross sections in NC DIS, measurements of the 
internal structure of jets in photoproduction and NC DIS and the study of the substructure 
dependence of jet cross sections in photoproduction. 

2 Parton evolution at low Bjorken-x 

To leading logarithm accuracy, the DGLAP evolution is equivalent to the exchange of a parton 
cascade with the exchanged partons strongly ordered in virtuality (Fig. |2J). The DGLAP 
equations sum the leading powers of a s log Q 2 in the region of strongly-ordered transverse 
momenta. However, DGLAP evolution is expected to breakdown at low x since only leading 
logarithms in Q 2 are resummed and contributions from log 1/x are neglected. These terms need 
to be taken into account since when logQ 2 <C log terms proportional to a s log 1/x become 
important. 




Figure 2: Example of Feynman diagram for forward jet production in NC DIS at low x. 

Several theoretical approaches exist which account for low-x effects not incorporated into the 
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DGLAP evolution. They are: (i) the BFKL evolution which resums large log(l/x) terms to all 
orders, this approach works at very low x and presents no kx ordering, the integration is taken 
over the full ft? phase scape of the gluons and the calculations make use of off-shell matrix 
elements together with unintegrated PDFs; (ii) the CCFM ^U] evolution provides angular- 
ordered parton emission and works for low and larger x, it is equivalent to BFKL for x — > 
and to DGLAP at large x\ (hi) the virtual-photon structure mimicks higher-order QCD effects 
at low x by introducing a second /^-ordered parton cascade on the photon side ^T] a la 
DGLAP, the resolved contribution is expected to contribute for (E^ 1 ) 2 > Q 2 and suppressed 
with increasing Q 2 . 

There exist several programs to calculate pQCD predictions of jet cross sections in NC DIS, eg 
Disent and Nlojet [T^j. These programs use the DGLAP evolution equations. Higher- 
order effects can be mimicked by the parton-shower approach in leading-logarithm Monte Carlo 
models like Rapgap ^3], which includes direct or direct plus resolved processes in virtual- 
photon interactions. Low-x effects not included in the DGLAP evolution are incorporated 
in Monte Carlo models like Cascade [T3| and Ariadne [16J. Cascade is based on the 
fcT-factorised unintegrated parton distributions and follows the CCFM evolution, whereas Ari- 
adne generates non-fc^ ordered parton cascades based on the color dipole model (CDM). 

Experimentally, deviations from the DGLAP evolution can be expected at low x and forward- 
jet rapidity since parton emission along the exchanged gluon ladder (see Fig. EJ) increases with 
decreasing x. Another method to obtain evidence for DGLAP breakdown is to study the 
azimuthal correlation between the two hardest jets. In DGLAP, partons entering the hard 
process with negligible kr produce a back-to-back configuration at LO. Values of A0 < 180° 
occur in DGLAP due to higher-order QCD effects. In models which predict a significant 
proportion of partons entering the hard process with large kx, the number of events with small 
A(p will increase. 

2.1 x dependence of the forward-jet cross section 

The forward-jet cross section has been measured for jets identified with the kx cluster 
algorithm in the longitudinally inclusive mode in the laboratory frame. Events with at least 
one jet of transverse energy -E^l AB > 3.5 GeV and 1.7 < rf^ B < 2.8 were selected. The 
events are required to fulfill the additional conditions: (i) Xj e t = Ej et /E p > 0.035, where E p is 
the proton-beam energy, and (ii) 0.5 < (E^ ) 2 /Q 2 < 5, following the proposal of Mueller and 
Navelet jTHj such as to allow for a large-level arm of evolution in x and to restrict the evolution 
in Q 2 . The measurements were made in the kinematic region given by 5 < Q 2 < 85 GeV 2 and 
0.0001 < x < 0.004. 

Figure El shows the forward-jet cross section as a function of x. The measured cross section rises 
with decreasing x. The NLO calculation corrected for hadronisation effects obtained using the 
program Disent with the renormalisation scale jj? R = 45 GeV 2 , the average (E 3 ^) 2 of the data, 
and the CTEQ6 parametrisations of the proton PDFs, is compared to the data in Fig. Et- The 
measured cross section is well described by the prediction for large values of x, but at low x 
values there is a large excess of data with respect to the calculation. This prediction is based 
on DGLAP evolution and so it is not expected to work in this region of phase space. 

Figure OJd shows the comparison with the predictions of different leading-logarithm parton- 
shower Monte Carlo models. The prediction of Rapgap including only direct processes is 
similar to the NLO calculation. Once the contribution from resolved processes is included, a 
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Figure 3: Measured forward-jet cross section (dots) [T7j as a function of Bjorken x. For compar- 
ison, the predictions of (a) Disent and (b) Rapgap (RG), Ariadne (CDM) and Cascade 
are also included. 



better description of the data is obtained for x down to 0.001. The Cascade prediction, based 
on the CCFM evolution, does not reproduce the shape of the data, whereas the prediction from 
CDM describes the data for x > 0.0015. In conclusion, no model can describe the sharp rise of 
the data at very low x, which remains a challenge and demands improved approximations to 
QCD in that region of phase space. 



2.2 Azimuthal jet separation 

The azimuthal correlation between the two hard jets in dijet events has been measured JH] 
using the fc^-cluster algorithm in the longitudinally inclusive mode in the 7*p centre-of-mass 
frame. The measurements were made in the kinematic region given by 5 < Q 2 < 100 GeV 2 
and 10~ 4 < x < 10 -2 . The cross sections refer to jets of > 5 GeV, — 1 < rf^ AB < 2.5 
and E^ max > 7 GeV, where E^ is the jet transverse energy in the 7*p centre-of-mass frame. 
Figured shows the measured distribution as a function of the azimuthal separation in different 
Q 2 regions. A significant fraction of events is observed at small azimuthal separations. Since a 
measurement of a multi-differential cross section as a function of x, Q 2 and A0* would be very 
difficult due to large migrations, the fraction of the number of dijet events with an azimuthal 
separation between and a, where a was taken as a = ~7r, was measured instead. The fraction 
S, defined as 

s _ f Q a N 2iet (A^,x,Q 2 )dAr 

/;iV 2jct (A0*,x,Q 2 )rfA0*' 

is better suited to test small-x effects than a triple-differential cross section. 

The measured fraction S as a function of x in different regions of Q 2 is presented in Fig. |3Jd. 
The data rise towards low x values, especially at low Q 2 . The predictions from Disent, which 
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Figure 4: (a) Measured A0* distribution in different regions of Q 2 (dots) For comparison, 
the predictions of Rapgap and Ariadne are also included, (b) Measured fraction S as a 
function of x in different Q 2 regions (dots) ^Hj- For comparison, the predictions of the fixed- 
order QCD calculations from Disent and Nlojet are also included. 



contain the lowest-order contribution to S, are several standard deviations below the data and 
show no dependence with x. On the other hand, the predictions of Nlojet, which incorporate 
NLO corrections to S, provide a good description of the data at large Q 2 and large x. However, 
they fail to describe the increase of the data towards low x values, especially at low Q 2 . This 
shows the need to incorporate NLO corrections to the calculations. 

Higher-order effects can be mimicked by the parton shower approach in Monte Carlo models. 
Figure shows the comparison with the predictions of Rapgap with direct only and resolved 
plus direct processes. A good description of the data is obtained at large Q 2 and large x. 
However, there is a failure to describe the strong rise of the data towards low x, especially at 
low Q 2 , even when including a possible contribution from resolved virtual-photon processes, 
though the description in other regions is improved. 

If the observed discrepancies are due to the influence of non-ordered parton emission, models 
based on the CDM or the CCFM evolution could provide a better description of the data. 
Figure |Sb shows the data compared with Ariadne and two predictions of Cascade, which 
use different sets of unintegrated parton distributions. These sets differ in the way the small-fc^ 
region is treated: in Jung2003 the full splitting function, i.e. including the non-singular term, 
is used in contrast to JS2001, for which only the singular term was considered. The predictions 
of Ariadne give a good description of the data at low x and Q 2 , but fail to describe the data 
at high Q 2 . The predictions of Cascade using JS2001 ^5] lie significantly above the data 
in all x and Q 2 regions, whereas those using Jung2003 [2U| are c l° ser to the data. Therefore, 
the measurement of the fraction S is sensitive to and can be used to gain information on the 
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Figure 5: Measured fraction S as a function of x in different Q 2 regions (dots) |19] . For com- 
parison, the predictions from (a) Rapgap and (b) Cascade and Ariadne are also included. 



unintegrated parton distributions. 



3 Multi-jet production in NC DIS 

Three-jet production in NC DIS provides a test of pQCD directly beyond LO since the since 
the lowest-order contribution is proportional to a 2 . Three-jet events arise from additional gluon 
brehmstrahlung or splitting of a gluon into a qq pair (Fig. EJ). 

The dijet (three-jet) cross section has been measured |2T] using the kp cluster algorithm in 
the longitudinally invariant mode in the Breit frame for events with at least two (three) jets 
of E^ B > 5 GeV and -1 < f?LAB < 2.5, and with dijet (three-jet) invariant masses in excess 
of 25 GeV; the kinematic region is defined by 10 < Q 2 < 5000 GeV 2 . Figure shows the 
dijet and three-jet cross sections as a function of Q 2 . The data are compared to the predictions 
of Nlojet up to 0(a 2 ) and 0(a 3 s ), respectively, using /i| = Q 2 + (E^ B ) 2 , where E^ B is 
the average jet transverse energy of the two (three) jets, and the factorisation scale \ip = Q. 
The CTEQ6 parametrisations of the proton PDFs have been used for the proton PDFs. The 
measured cross sections are well described by the predictions. 

The Q 2 dependence of the ratio of the three-jet to the dijet cross section has been studied (see 
Fig. |7|d). Many experimental and theoretical uncertainties cancel in the ratio and therefore 
this observable provides a more accurate test of color dynamics, especially at low Q 2 since the 
theoretical uncertainty of the ratio is of the same order as at higher Q 2 . The calculations from 
Nlojet, using the five sets of the CTEQ4 "A-series" of proton PDFs, give a good description 
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Figure 6: Examples of Feynman diagrams for three-jet production in NC DIS: (a) gluon brehm- 
strahlung from a quark and (b) gluon splitting into a qq pair. 

of the data. The predictions for different values of a s show the sensitivity of this observable to 
the value of a s (Mz) assumed in the calculations. A value of 

a s (M z ) = 0.1179 ±0.0013 (stat.) +°;°046 (exp.) t° f t\ (th.) 

has been extracted from the ratio, which is in good agreement with the world average, 0.1182 ± 
0.0027 [22] and other determinations of a s from HERA data. 



4 Jet substructure and the dynamics of quarks and glu- 
ons 

The internal structure of a jet depends mainly on the type of primary parton -quark or gluon- 
from which it originated and to a lesser extent on the particular hard scattering process. QCD 
predicts that at sufficiently high , where fragmentation effects become negligible, the jet 
structure is driven by gluon emission off the primary parton and is then calculable in pQCD. 
QCD also predicts that gluon jets are broader than quark jets due to the larger colour charge 
of the gluon. This prediction has been confirmed at LEP in measurements of the internal 
structure of jets [22]. 

The internal structure of jets has been studied by means of the jet shape. The integrated jet 
shape is defined as the fraction of the jet transverse energy that lies inside a cone in the rj — (p 
plane of radius r concentric with the jet axis, using only those particles belonging to the jet [24J: 

where Et{t) is the transverse energy within the given cone of radius r. The mean integrated 
jet shape, (ip(r)), is defined as the averaged fraction of the jet transverse energy inside the cone 
of radius r: 

E T (r) 



s 



ZEUS ZEUS 




Figure 7: (a) Measured dijet and three-jet cross sections in NC DIS as a function of Q 2 
(dots) [21] • (b) Measured ratio of the three-jet to the dijet cross section as a function of 
Q 2 (dots) [21]. For comparison, the predictions of Nlojet are also included in both figures. 



where iVj ets is the total number of jets in the sample. 

The jets have been identified using the kx cluster algorithm in the longitudinally inclusive mode 
and selected according to E 3 ^ 1 > 17 GeV and —1 < rj^ ct < 2.5. The kinematic region in the 
photoproduction sample is defined by Q 2 < 1 GeV 2 and 142 < W 1P < 293 GeV, where W 1P is 
the 7p centre-of-mass energy, and in the NC DIS sample by Q 2 > 125 GeV 2 . 

4.1 Jet-shape measurements 

The measured mean integrated jet shape as a function of r, (ip(r)), for different regions in rf et 
is shown in Fig. [H^ for the photoproduction regime [23]. The jets broaden as rf et increases. 
Leading-logarithm parton-shower predictions from Pythia [20] for resolved plus direct pro- 
cesses and for gluon- and quark-initiated jets are compared to the data in Fig. [H^l. The mea- 
sured (ip{r)) is reasonably well described by the MC calculations of Pythia for resolved and 
direct processes for —1 < ?7 Jet < 1.5, whereas for 1.5 < rf ct < 2.5, the measured jets are slightly 
broader than the predictions. From the comparison with the predictions for gluon- and quark- 
initiated jets, it is seen that the measured jets are quark-like for —1 < rf et < and become 
increasingly more gluon-like as if et increases. 

The measured (ip{r)} for different regions of Ej, is shown in Fig. [Sjo for NC DIS events. The 
NLO QCD calculations of (ip(r)), corrected for hadronisation and Z°-exchange effects, are 
compared to the data in the figure. The NLO QCD calculations give a good description of 
(ip{r)) for r > 0.2; the fractional differences between the measurements and the predictions are 
less than 0.2% for r = 0.5. 

The quark and gluon content of the final state has been investigated in more detail by studying 



9 



ZEUS 



ZEUS 



i r 

0.8 r 

0.6 7 

0.4 ~ T 

0.2 ~ 7 

- 





1 < if* < 

•ZEUS 98-00 TP i ,j Hadron level 
PYTHIA 
gluons 

quarks 0.1 f 

_Q Q o O O O O O - 2 Q Q O o O O 2 



o (data-PYTHIA)/PYTHIA- 




- 5 o o o o 
0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



v 0.5 



v 0.5 



v 0.5 





u 

v 0.5 



1 | , 1 J , 1 , 1 1 | 1 J I III! 

-I< rf* < 2.5^_^-»— •— ■— ■ ~ 

17 < E£'< 21 GeV 

~ J ~- 
- m 

■ ■ ZEUS 98-00 DIS 


9^ 21 < E£'< 25 GeV 
7 _/ Hadron level 
: □ (data-DISENT)/DLSENT : 




" 1 1 1 1 1 1 1 1 1 1 II 1 1 1 1 1+ 


" i i i i i i i i i i i ihiji-i- 




I f 25 < E*'< 29 GeV \ 

- ■ — DISENT 

- n 


I 29 < E*'< 35 GeV I 

- □ Q n q - 


1 1 1 1 1 1 1 1 II ill 1 h- 


1 1 i i i i 1 1 ii ih i h- 


/ 35 < E*'< 41 GeV I 


I 41 < E^r'< 47 GeV I 


'l 1 1 1 1 1 1 1 1 1 1 ill 1 l-H 


" 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 l-H 


: ./ 47 < E'«< 55 GeV I 


I S 55 < E'«< 71 GeV I 







0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

r r 



(a) (b) 

Figure 8: (a) Measured mean integrated jet shape as a function of r in different regions of r)^ et 
in photoproduction (dots) |25] . (b) Measured mean integrated jet shape as a function of r in 
different regions in NC DIS (squares) For comparison, the predictions of (a) Pythia 
and (b) Disent are also included. 



the ?7j et dependence of the mean integrated jet shape in photoproduction and NC DIS at a 
fixed value of r = 0.5 (see Fig. Et)- The jet shape at a fixed value of r = 0.5 decreases 
with increasing 77^* for photoproduction, whereas the jets in NC DIS show no dependence with 
rf et . The comparison of the data with the predictions for quark and gluon jets shows that the 
NC DIS jets are consistent with being dominated by quark jets, whereas the broadening of 
the jets in photoproduction is consistent with an increase of the fraction of gluon jets as rf et 
increases. The dependence of the mean integrated jet shape for a fixed value of r = 0.5 as a 
function (see FigEb) in NC DIS shows that the jets become narrower as E 3 ^ increases. 
The same is observed in photoproduction (not shown). The comparison of the data with the 
NLO calculations assuming different values of a s (Mz) shows the sensitivity of this observable 
to the value of a s (M z ). A value of a s of 

a s (M z ) = 0.1176 ±0.0009 (stat.) j™™ (exp.) t° Zl (th.) 

was determined from this observable. This determination of a s has experimental uncertainties 
as small as those based on previous measurements. However, the theoretical uncertainty is large 
and dominated by terms beyond NLO. Further theoretical work on higher-order contributions 
would allow an improved determination of a s from the integrated jet shape in DIS. 



4.2 Substructure dependence of jet cross sections 

The predictions of the Monte Carlo for the jet shape reproduce well the data and show the 
expected differences for quark- and gluon-initiated jets. These differences are used to select 
samples enriched in quark and gluon jets to study in more detail the dynamics of the hard 
subprocesses. The predicted shapes of the distribution in ip(r = 0.3) for quark and gluon jets 
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Figure 9: (a) Measured mean integrated jet shape at a fixed value of r = 0.5 as a function of 
rf et in photoproduction (dots) and NC DIS (squares) [22] • For comparison, the Monte Carlo 
predictions for quark- and gluon-initiated jets are also included, (b) Measured mean integrated 
jet shape at a fixed value of r = 0.5 function of £$* in NC DIS fZ%. For comparison, the 
predictions of Disent assuming three different values of a s are also included, (c) Predicted 
integrated jet shape distributions from Pythia at r = 0.3 for quark- (hatched histogram) and 
gluon-initiated (shaded histogram) jets. 

are different, as shown in Fig. Eh- A sample enriched in quark (or "narrow") jets was selected 
by requiring an integrated jet shape above 0.8 and a sample enriched in gluon (or "broad") jets 
was selected by requiring an integrated jet shape below 0.6. Pythia predicts a purity of 57% 
for gluons and 84% for quarks, and the efficiencies are 58% for gluons and 51% for quarks. 

The differential inclusive-jet cross-section da/drf et for photoproduction is shown in Fig. fTUk 
for samples of broad and narrow jets, separated according to the selection explained above. 
The measured cross sections exhibit different behaviours: the rf et distribution for broad jets 
increases up to the highest rf^ value measured, whereas the distribution for narrow jets peaks 
at rf ct « 0.7. Monte Carlo calculations using Pythia are compared to the measurements in 
Fig.llOa. The same selection method was applied to the jets of hadrons in the MC event samples 
and the calculations have been normalised to the total measured cross section of each sample. 
The MC predictions provide a good description of the shape of the narrow-jet distribution in 
the data. The shape of the broad-jet distribution in the data is reasonably well described by 
Pythia. From the calculation of Pythia, the sample of broad jets selected according to the 
jet shape is predicted to contain 15% of gg subprocesses in the final state and 50% of gq, and 
a contamination from processes with only quarks in the final state of 35%. There is a large 
contribution from gq final states in the broad-jet sample because the partonic cross section for 
the resolved subprocess q^g v — > qg is much larger than the cross section for the subprocesses 
QQ — > 99 P ms 99 ~ ¥ 99- The sample of narrow jets contains 62% of qq subprocesses and 
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Figure 10: (a) Measured inclusive-jet cross section in photoproduction as a function of rf et for 
samples of broad (dots) and narrow (open circles) jets j2Hj. (b) Measured dijet cross section 
in photoproduction as a function of cos#£ road (stars) |2Sj- For comparison, the predictions of 
Pythia for jets selected in the same way as in the data and for quark- and gluon-initiated jets 
are included. 



34% of qg, with a contamination of 4% from gg subprocesses. Figure ITUk also shows the 
predictions of Pythia for jets of quarks and gluons separately. These predictions have been 
obtained without any selection and are normalised to the data cross sections. The calculation 
that includes only quark-initiated jets gives a good description of the narrow-jet cross section, 
whereas the calculation for gluon-initiated jets provides a reasonable description of the broad- 
jet cross section. This result supports the expectation that the broad (narrow)-jet sample is 
dominated by gluon (quark)-initiated jets. 

The distribution in 9*, where 9* is the angle between the jet-jet axis and the beam direction 
in the dijet system, reflects the underlying parton dynamics and is sensitive to the spin of the 
exchanged particle. In the case of direct-photon interactions, the contributing subprocesses 
at LO QCD involve quark exchange and so da/d | cos 9*\ oc (1 — |cos#*|) _1 as |cos#*| — > 1. 
In the case of resolved-photon interactions, the dominant subprocesses are those that involve 
gluon exchange and da/d\ cos6>*| oc (1 — | cos9*\)~ 2 as | cos#*| — > 1. The study of the angular 
distribution for dijet events with tagged quark- and/or gluon-initiated jets in the final state, 
provides then a handle to investigate the underlying parton dynamics further. 

The sample of photoproduced dijet events with one broad jet and one narrow jet was used to 
measure da/d cos 9^ Toad , where 0£ road re f ers to the scattering angle measured with respect to 
the broad jet. Figure UTIb shows the measured dijet cross section as a function of cos6^ road . 
The measured and predicted cross sections were normalised to unity at cos#£ road = 0.1. The 
dijet angular distribution shows a different behaviour on the negative and positive sides; the 
measured cross section at cos#£ road = 0.7 is approximately twice as large as at cos#£ road = 
—0.7. The calculation from Pythia gives a good description of the shape of the measured 
da/d cos #broad- The predictions of Pythia for the partonic content are: 52% of qg subprocesses, 
4% of gg and 44% of qq. The observed asymmetry is adequately reproduced by the calculation 
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and is understood in terms of the dominant resolved subprocess q^g p — > qg. The cos0£ road 
distribution for this subprocess is asymmetric due to the different dominant diagrams in the 
regions cos0£ road — > ±1: t-channel gluon exchange (cos#£ road — > +1) and w-channel quark 
exchange (cos^ road -> -1). 

In conclusion, the hard subprocesses have been investigated separately in photoproduction for 
the first time using the internal structure of jets. 



5 Conclusions 

HERA has become a unique QCD-testing machine. At large scales considerable progress in 
understanding and reducing the experimental and theoretical uncertainties has led to very 
precise measurements of the fundamental parameter of the theory, the strong coupling constant 
a s (see Fig. ITTj) as well as further insight into the dynamics of quarks and gluons. The use 
of observables based on the application of jet algorithms to the hadronic final state of deep 
inelastic scattering and of photon-proton interactions leads now to determinations that are as 
precise as those coming from more inclusive measurements, such as from r decays. To obtain 
even better accuracy in the determination of a s , further improvements in the QCD calculations 
are needed, e.g. next-to-next-to-leading-order corrections. 



th. uncert. 



exp. uncert. 



Jet shapes in NC DIS 

ZEUS (DESY 04-072 - hep-ex/0405065) 

Multi-jets in NC DIS 

ZEUS prel. (contributed paper to ICHEP04) 
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ZEUS (Phys Lett B 560 (2003) 7) 

Subjet multiplicity in CC DIS 

ZEUS (Eur Phys Jour C 31 (2003) 149) 

Subjet multiplicity in NC DIS 

ZEUS (Phys Lett B 558 (2003) 41) 

NLO QCD fit 

HI (Eur Phys J C 21 (2001) 33) 
NLO QCD fit 

ZEUS (Phys Rev D 67 (2003) 012007) 

Inclusive jet cross sections in NC DIS 
HI (Eur Phys J C 19 (2001) 289) 
Inclusive jet cross sections in NC DIS 
ZEUS (Phys Lett B 547 (2002) 164) 

Dijet cross sections in NC DIS 
ZEUS (Phys Lett B 507 (2001) 70) 

World average 

(S. Bethke, hep-ex/0407021) 



0.1 0.12 



0.14 

a s (M z ) 



Figure 11: Summary of a s determinations at HERA compared with the world average. 



At low x considerable progress has also been obtained in understanding the mechanisms of par- 
ton emission, though the interplay between the DGLAP, BFKL and CCFM evolution schemes 
has still to be fully worked out. Further progress in this respect needs both more experimental 
and more theoretical work. 
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